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[ Context. Classical Be stars are hot non-supergiant stars surrounded by a gaseous circumstellar disk that is responsible for the observed 

[ I , infrared-excess and emission lines. The phenomena involved in the disk formation still remain highly debated. 

Aims. To progress in the understanding of the physical process or processes responsible for the mass ejections and test the hypothesis 
that they depend on the stellar parameters, we initiated a survey on the circumstellar environment of the brightest Be stars. 
Methods. To achieve this goal, we used spectro-interferometry, the only technique that combines high spectral (R= 12000) and high 
spatial (0 min =4mas) resolutions. Observations were carried out at the Paranal observatory with the VLTI/ AMBER instrument. We 
concentrated our observations on the Bry emission line to be able to study the kinematics within the circumstellar disk. Our sample is 
composed of eight bright classical Be stars : a Col, k CMa, oj Car, p Car, 6 Cen, \i Cen, a Ara, and o Aqr. 

Results. We managed to determine the disk extension in the line and the nearby continuum for most targets. We also constrained the 
disk kinematics, showing that it is dominated by rotation with a rotation law close to the Keplerian one. Our survey also suggests that 
I ■ these stars are rotating at a mean velocity of V/V c = 0.82 ± 0.08. This corresponds to a rotational rate of £2/f2 c = 0.95 ± 0.02 

Conclusions. We did not detect any correlation between the stellar parameters and the structure of the circumstellar environment. 
Moreover, it seems that a simple model of a geometrically thin Keplerian disk can explain most of our spectrally resolved K-band 
data. Nevertheless, some small departures from this model have been detected for at least two objects (i.e, k CMa and a Col). Finally, 
our Be stars sample suggests that rotation is the main physical process driving the mass-ejection. Nevertheless, smaller effects from 
other mechanisms have to be taken into account to fully explain how the residual gravity is compensated. 

Key words. Techniques: high angular resolution - Techniques: interferometric - Stars: emission-line, Be - Stars: winds, outflows - 
Stars: individual (a Col, k CMa, a> Car, p Car, S Cen, ji Cen, a Ara, o Aqr) - Stars: circumstellar matter 
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. 1 . Introduction large spread of rotational velocities, whereas the cooler Be stars 

are more likely to be critical rotators. Moreover, Abbott (1979) 
Classical Be stars are hot non-supergiant stars that have at least showed that the radiative pressure alone can initiate mass ejec- 
once exhibited the so called "Be-phenomenon", i.e. Balmer lines tion only for the earliest Be stars Consequently, it is not yet clear 
in emission and IR-excess originating from a dense gaseous if Be stars can be consider ed as a homogeneous group of stars in 
circumstellar environment. A generally accepted scheme is the terms of mass . e j e ction processes (Stee & Meilland 2009). 
presence of two distinct regions in the envelope: a dense equa- 
torial disk dominated by rotation that is responsible for most of An efficient way to test these hypotheses is to constrain the 
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the line emission and IR-excess, and a more diluted radiatively geometry and kinematics of the Be stars circumstellar environ- 
ed ! driven polar wind (Lamers & Waters 1987) with terminal veloc- ments because they mainly depend on the mechanisms respon- 
■ - ' ities on the order of several hundreds of km.s-' (Marlborough sible for the mass-ejection. This can only be performed with 
1982). long-baseline interferometry with sufficient spectral resolution, 
However, the physical process or processes responsible for as shown b ? Meilland et al. (2007a ; 2007b; 2011) using the 
the mass ejection and reorganization of matter in the circumstel- VLTI/ AMBER on the stars a Ara, k CMa, and 5 Sco, or Delaa et 
lar environment are still highly debated. The relative effect of ro- aL ( 201 °) usln 8 the CHARA/VEGA instrument on 48 Per and 0> 
tation, radiative pressure, pulsation, and binarity still need to be Per In a few cases ' the stellar photosphere can also be resolved 
quantified. For instance, classical Be stars are known to be fast- b y the interferometer, allowing measuring the photosphere's flat- 
rotators, but estimates of their rotational velocities range from tenm g- In thls case the physical stellar parameters can be deter- 
50 to 100% of the critical velocity which means that rotation mmed ' as for Achernar (Domiciano de Souza et al. 2003). For 
alone might not explain the mass-ejection for all cases. Cranmer th i s specific ob j ect < a clear signature of a polar wind was also 
(2005) showed in a statistical study of the velocity of 462 classi- detected using the same VLTI/VINCI dataset by Kervella et al. 
cal Be stars that the hottest Be stars (i.e., T eff >18000K) have a (2006). 

Before the availability of a new generation of instruments 
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Table 1. Be stars observed in this survey and their stellar parameters taken from the literature or estimated from the SED fit. 



Name 


HD 


m K 


spectral 


distance 


T e ff 


vsini 


v c 


pol. angle 


estim. 


estim. F t ,„,,, 






(mag) 


class 


(pc) 


(K) 


(krns" 1 ) 


(kms -1 ) 


(dog) 


(Re) 


(in the K band) 


a Col 


37795 


2.8 


B7IV 


80±2 


12963+203 


192+12 


355+23 


109 


5.8 


0.25 


KCMa 


50013 


3.5 


B1.5V 


202±5 


24627+590 


244+17 


535+39 


106 


5.9 


0.47 


oj Car 


89080 


3.5 


B8III 


105+1 


11720+431 


245+13 


320+17 


38 


6.2 


0.20 


p Car 


91465 


3.2 


B4V 


148±9 


17389+415 


285+20 


401+28 


68 


6.0 


0.45 


(5 Cen 


105435 


2.7 


B2IV 


127+8 


22360+589 


263+14 


527+29 


137 


6.5 


0.45 


H Cen 


120324 


4.0 


B2IV-V 


155+4 


22554+661 


166+10 


508+32 




5.5 


0.37 


a Ara 


158427 


2.5 


B3IV 


82±5 


18044+310 


305+15 


477+24 


166 


5.5 


0.56 


o Aqr 


209409 


4.6 


B7IV 


133+4 


12942+402 


282+20 


391+27 


6 


4.0 


0.31 



studies were conducted in photometric bands to measure exten- 
sion and flattening of the circumtellar environments of many 
Be stars. For instance, y Cas and if/ Per were observed by 



Tycner et al. (2006) with the NPOI interferometer in the Ha do- 



main. They found that a uniform disk or a ring-like model were 
inconsistent with their data, whereas a Gaussian model well fit- 
ted the measurements, y Cas' disk was also consistent with the 
orbital parameters already published. Nevertheless, higher preci- 
sion binary solutions were mandatory to test for a possible disk 
truncation by the secondary. The disk of if/ Per was found to be 
truncated by a companion as already predicted by Waters ( 1 986) 
y Cas, ijj Per, £ Tau and k Dra were also observed with 
the CHARA interferometer in the K-band by |Gies etal. (2007)[ 
Using Gaussian elliptical fits of visibilities, they found that the 
disk size in the K-band was smaller than in Ha owing to a higher 
Ha opacity and relatively higher neutral hydrogen fraction with 
increasing disk radius. All these Be stars are known binaries, and 
this binarity effect was found to be more significant for i// Per and 
k Dra. 



Tycner et al. (2008) observed^ Oph with the NPOI interfer- 



ometer and obtained a good fit of the Hot emitting disk with a 
circularly symmetric Gaussian, favoring the hypothesis that this 
object is seen under a very low inclination angle. 

Using the mid-infrared VLTI/MIDI instrument, 
Meill and et al. (2009)| have observed seven classical Be 
stars between 8 and 12/rni. They found that the size of the disk 
does not vary strongly with wavelength within this spectral 
domain, which is a very different conclusion compared to the 
results for B[e] stars where increasing sizes as a function of 
wavelength where found (Millour et al 2009; Meilland et al. 
2010; Borges Fernandes et al. 2011). Moreover, the size of 
a Arae's disk was found to be identical at 2, 8, and 12 fim, 
which might be because of to a disk truncation by a companion. 
Finally, it seems from Meilland et al.'s (2009) studies that the 
disks of late type Be stars might be smaller than those of the 
early type. 

Long-baseline interferometry is also a powerful technique 
for detecting companions, which was shown in Meilland et al. 
(2008) for the Be star 6 Cen and Millour et al. (2009) for the B[e] 
star HD 87643. Moreover, in a theoretical study of the forma- 
tion and dissipation of the Be stars' equatorial disks, Meilland et 
al. (2006) showed that interferometric follow-up of these events 
is the best-suited technique to deduce the physical parameters 
of the system. However, Kanaan et al. (2008) showed that, for 
Achernar, coupling spectroscopic follow-up and large-band in- 
terferometric observations at one epoch was sufficient to roughly 
understand the geometry and kinematics of this star. 

To progress in the understanding of Be stars, we initiated 
an observational campaign of the brightest, closest objects us- 
ing the VLTI/AMBER (Petrov et al. 2007; Robbe-Dubois et al. 
2007) and VLTI/MIDI (Leinert et al. 2003) instruments for the 



southern stars and the CHARA/VEGA (Mourard et al. 2009) for 
the northern ones. In this paper, we present new VLTI/AMBER 
spectro-interferometric observations of eight classical Be stars : 
a Col, k CMa, a> Car, p Car, 6 Cen, fi Cen, a Ara, and o Aqr. 

The paper is organized as follows. In Sect. 2 we briefly in- 
troduce each target and constrain their physical parameters from 
various sources in the literature. The observations and data re- 
duction process are then presented in Sect. 3, and in Sect. 4 we 
provide a qualitative analysis of the reduced data for each object. 
Then, in Sect. 5, the data are analyzed with various "toy mod- 
els". Finally, these results are discussed in Sect. 6 and a short 
conclusion is drawn in Sect. 7. 



2. Our Be star sample 

Considering the actual limiting magnitude of the VLTI/AMBER 
instrument, i.e. H=K=5 for an unresolved source observed with 
the 1.8m auxiliary telescopes in medium or high spectral resolu- 
tion modes, about 30 classical Be stars are observable. However, 
to constrain the circumstellar disk kinematics efficiently, we de- 
cided to exclude targets with weak emission lines or transient 
disks. Moreover, for this first survey we mainly concentrated on 
targets with mx <4. Finally, we also decided to include our first 
attempt to observed a fainter target, i.e. o Aqr (m^ =4.6). 

Because one of our goals is to determine whether or not 
the Be phenomenon depends on the basic stellar parameters, we 
tried to select targets with the widest range of spectral classes 
possible. Finally, the stars selected in our sample, i.e. a Col, k 
CMa, lo Car, p Car, 5 Cen, fi Cen, a Ara, and o Aqr, have spec- 
tral types ranging from B1.5 to B8, and luminosity class of III, 
IV and V. TableQ]briefly describes the targets: 

• Spectral class and the K-band magnitude (mK) are taken 
from the CDS databases. 

• Distance (d) is derived from van Leeuwen (2007). 

• Effective temperature (T e ff), v sini, and critical velocity (V c ) 
are taken from Fremat et al. (2005) 

• Polarization measurement is taken from Yudin (2001) 

• Stellar radius (R*) and K-band environment relative flux 
(F env ) are estimated by fitting the spectral energy distribu- 
tion (SED) using reddened Kurucz (1979) models for stel- 
lar atmospheres and using stellar parameters (T e ff and g e ff) 
from Fremat et al. (2005). The SED is first reconstructed 
using photometric and spectro-photometric measurements 
from the ultraviolet (IUE specta) to the far-infrared (IRAS 
data). To avoid contamination from the circumstellar flux, 
the fit of the stellar contribution to the flux is made from the 
ultraviolet to the visible. (see Meilland et al. 2009 for more 
details). 
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Table 2. VLTI/ AMBER observing log. 



Obs. time 




Telescopes 


Length 




Position angle 


Instrument 


DIT 


Coherence 


Seeing 


Calibrators 


(UTC) 




conf. 


(m) 




O 


mode 


(s) 


(ms) 


o 


(HD) 










a 


Col 












2008-01-06 03:39 


K0-G1-A0 


88.6/90.5/126.2 




-150.1/ -59.7/- 104.3 


LR-K-F 


0.05 


5.4 


1.10 


81188 


2010-01-09 00:48 


D0-H0-K0 


60.2/30.1/ 90.3 




51.4/51.4/51.4 


HR-K-F 


5.00 


3.1 


1.25 


34642 


2010-01-12 00:48 


G1-D0-H0 


27.3/ 50.0/ 66.6 




-68.0/ -3.8/ -25.4 


HR-K-F 


5.00 


3.9 


1.06 




2010-01-20 01:35 


K0-G1-A0 


90.0/ 88.7/127.5 




-157.9/ -68.9/-1 13.8 


HR-K-F 


5.00 


5.3 


0.85 


34642 


2010-01-20 01:56 


K0-G1-A0 


89.9/ 89.6/128.0 




-156.1/ -67. 1/-111.7 


HR-K-F 


5.00 


3.8 


1.30 


34642 


kCMsl 


2008-12-18 03:34 


U1-U3-U4 


102.2/ 55.6/125.3 




17.1/ 96.1/ 42.9 


HR-K-F 


1.00 


5.3 


0.63 




2008-12-18 07:32 


U1-U3-U4 


94.5/58.9/115.3 




43.9/ 129.2/ 74.4 


HR-K-F 


1.00 


4.9 


0.64 




2008-12-20 03:37 


K0-G1-A0 


90.5/81.4/121.3 




-167.9/ -77.5/-125.7 


LR-K-F 


0.25 


4.4 


0.78 


40805 


2008-12-24 03:50 


K0-G1-A0 


90.5/ 84.8/124.2 




-164.3/ -74.5/-121.2 


LR-K-F 


0.05 


4.6 


0.79 


27442, 57299 


2008-12-24 05:07 


K0-G1-A0 


89.9/90.1/128.0 




-154.7/ -65.3/- 109.9 


LR-K-F 


0.05 


3.5 


1.00 


48305, 57299 


2008-12-24 08:01 


K0-G1-A0 


80.9/ 82.8/102.2 




-139.8/ -37.0/ -87.6 


LR-K-F 


0.05 


5.4 


0.66 


48305, 57299 


2010-01-09 01:39 


D0-H0-K0 


57.8/ 28.9/ 86.6 




46.2/ 46.2/ 46.2 


HR-K-F 


5.00 


3.1 


1.26 


34642 


2010-01-18 00:24 


H0-G0-E0 


27.4/ 13.7/41.0 




-142.8/-142.8/-142.8 


HR-K-F 


5.00 


3.6 


0.71 


54173 










CO 


Car 












2008-12-21 04 


49 


K0-G1-A0 


74.2/72.9/128.0 




175.8/- 125.0/- 154.9 


HR-K-F 


5.00 


3.8 


0.78 


75063, 69596 


2008-12-21 07 


37 


K0-G1-A0 


72.1/ 83.0/126.6 




-156.4/ -85.7/-1 18.2 


HR-K-F 


5.00 


3.4 


0.89 


69596 


2008-12-24 06 


33 


K0-G1-A0 


73.7/79.3/127.7 




-167.4/-100.5/-132.6 


LR-K-F 


0.05 


3.0 


1.20 


98134, 57299 










P 


Car 












2008-12-23 07:31 


K0-G1-A0 


79.3/82.3/127.8 




-161.2/ -85.7/-122.7 


MR-K-F 


1.00 


3.3 


1.11 


94286, 69596 


2009-03-22 03:59 


K0-G1-A0 


71.7/89.6/120.8 




-138.0/ -54.4/ -90.6 


MR-K-F 


1.00 


6.0 


0.70 


94286 


tfCen 


2009-03-21 04 


18 


K0-G1-A0 


83.9/87.9/127.4 




-153.8/ -69.5/- 110.4 


MR-K-F 


1.00 


10.9 


0.74 


110458 


2010-01-20 08 


12 


K0-G1-A0 


84.1/ 87.7/127.5 




-154.5/ -70.4/-1 11.4 


LR-HK 


0.05 


2.6 


1.27 


110458 


2011-05-19 02 


41 


U1-U2-U4 


46.4/82.3/111.4 




43.8/ 106.6/ 84.8 


LR-HK 


0.025 


1.3 


1.47 


103513 










/' 


Cen 












2011-06-26 23:21 


K0-A1-G1 


128.9/78.1/88.2 




-121.5/ 100.6/-157.9 


HR-K-F 


6.00 


1.6 


0.90 


128488 


2011-06-30 00:19 


D0-I1-G1 


82.3/ 44.5/ 70.8 




104.8/-134.5/ 137.5 


HR-K-F 


5.00 


3.4 


0.94 


128488 










a 


Ara 












2007-07-28 05 


51 


Gl-DO-HO 


71.5/ 44.7/55.4 




-5.4/ 123.9/33.3 


LR-K 


0.025 


3.5 


0.44 


177716, 164371 


2007-04-13 05 


47 


HO-GO-EO 


30.8/ 15.4/46.3 




-143.4/-143.4/-143.4 


LR-K 


0.025 


3.0 


0.63 


124454 


2007-06-09 07 


11 


K0-G1-A0 


74.2/90.0/112.3 




-135.8/ -41.6/ -82.8 


LR-K 


0.025 


2.7 


0.57 


166460, 164371 


2007-06-09 08 


02 


K0-G1-A0 


68.3/ 88.9/102.5 




-130.4/ -30.7/ -71.7 


LR-K 


0.025 


1.6 


0.95 


166460, 164371 


2007-06-06 07 


30 


K0-G1-A0 


72.5/ 89.6/109.5 




-134.2/ -38.5/ -79.7 


LR-K 


0.025 


1.3 


1.38 


124454, 164371 


2007-06-06 08 


50 


K0-G1-A0 


63.4/88.0/95.0 




-126.9/ -22.7/ -63.0 


LR-K 


0.025 


1.3 


1.42 


164371 


2007-04-14 06 


22 


HO-GO-EO 


31.4/ 15.7/47.2 




-133.5/-133.5/-133.5 


LR-K 


0.025 


2.9 


0.67 


164371 


2007-05-16 07 


34 


HO-DO-AO 


61.1/ 30.6/91.7 




-96.4/ -96.4/ -96.4 


LR-K 


0.025 


3.0 


1.25 


124454, 164371 


2007-06-10 07 


44 


K0-G1-A0 


70.2/ 89.2/105.7 




-132.0/ -34.1/ -75.2 


LR-K 


0.025 


2.3 


0.58 


164371, 21201 


2007-05-17 02 


34 


HO-DO-AO 


60.2/30.1/ 90.3 




-156. 3/- 156. 3/- 156. 3 


LR-K 


0.025 


3.0 


0.79 


124454 


2007-05-17 09 


53 


HO-DO-AO 


49.3/24.6/73.9 




-67.2/ -67.2/ -67.2 


LR-K 


0.025 


1.9 


1.53 


164371 


2011-06-30 02 


04 


D0-I1-G1 


78.2/ 45.8/ 62.7 




82.6/-150.7/ 118.4 


HR-K-F 


5.00 


3.7 


0.85 


163145, 152786 


2011-06-30 03 


03 


D0-I1-H0 


81.4/37.2/63.8 




93.9/ -36.5/ 67.6 


HR-K-F 


5.00 


2.9 


1.08 


163145, 152786 


2011-06-30 03 


37 


D0-I1-H0 


82.2/38.1/ 63.2 




100.4/ -31.9/ 73.9 


HR-K-F 


5.00 


2.9 


1.05 


163145, 152786 


2011-06-30 04 


28 


D0-I1-H0 


82.0/39.0/61.2 




110.1/ -24.7/83.3 


HR-K-F 


5.00 


2.4 


1.27 


163145, 152786 


2011-06-30 05 


03 


D0-I1-H0 


81.0/39.4/58.9 




117.5/ -19.2/ 90.2 


HR-K-F 


5.00 


2.9 


1.03 


163145, 152786 


2011-06-30 06 


00 


D0-I1-H0 


78.5/ 39.8/ 54.4 




129.4/ -10.7/ 101.3 


HR-K-F 


5.00 


2.9 


1.02 


163145, 152786 


2011-06-30 06 


38 


D0-I1-H0 


76.4/ 39.9/ 50.6 




138.1/ -4.7/ 109.7 


HR-K-F 


5.00 


3.2 


0.92 


163145, 152786 










o 


Aqr 












2011-06-30 10:08 


D0-I1-H0 


67.4/ 34.3/ 58.8 




102.8/ -16.5/ 72.2 


HR-K-F 


5.00 


2.9 


1.02 


209926 



3. Observations and data reduction process 

The VLTI/ AMBER observations of classical Be stars were car- 
ried out at Paranal Observatory between 2007 and 2011. Most 
observations were executed using the 1.8 meters movable aux- 
iliary telescopes (AT), a few others were made using 8.0m unit 
telescopes (UT). The VLTI/ AMBER observing log for the eight 
targets and their corresponding interferometric calibrators are 
presented in Table [2] The (u,v) plan coverage for each target is 
plotted in Fig[TJ 

Because the observations were spread over a long period 
of time, starting soon after the opening of AMBER to the sci- 



entific community, the data are quite inhomogeneous, both in 
terms of quality and observing modes. AMBER offers three dif- 
ferent spectral dispersion modes: the low-resolution (LR) with 
R=A/(M^35, the medium-resolution (MR) with R^1500, and the 
high-resolution (HR) with R^ 12000 . All these modes were used 
during our observing campaign, depending on their availability 
at the epoch of observations, on the brightness of the target, and 
on the seeing conditions. The MR and HR observations were 
centered on the Bry emission line to enable the study of the cir- 
cumstellar gas kinematics through the Doppler effect. 

Most of the data have benefited from the installation of 
the fringe tracker FINITO that enables longer integration time 
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Fig. 1. (u,v) plan coverage obtained for the observed Be stars. LR mode observations are plotted as squares, MR ones as triangles, 
and HR ones as circles. Each observation, i.e. three baseline measurements, is plotted with a different color. 



by stabilizing the fringes. The detector integration time (DIT) 
ranges from 25 ms for bright targets in LR mode without 
FINITO to several seconds for observations in HR mode with 
FINITO. Under good seeing conditions (i.e. seeing <0.8") 
FINITO allows one to obtain a significantly higher signal-to- 
noise ratio (S/N) both for short and long integration time, signifi- 
cantly enhancing the quality of the data for all observing modes. 

Data were reduced using the VLTI/ AMBER data reduction 
software, i.e., amdlib v3.8.3bl (see Tatulli et al. 2007 and 
Chelli et al. 2009 for detailed informations on the AMBER 
data reduction). We selected individual exposures with the stan- 
dard selection criteria (Millour et al. 2007). We rejected 80% of 
frames with the lowest S/N. For observations in LR mode we 
also rejected the frames with a piston larger than lOyimi as well 
as frames with a flux ratio between the beams higher than 3. 

The interferometric observables (visibility, differential 
phase, and closure phase) were then averaged and calibrated. For 
this last step we used scripts described in Millour et al. (2007) 
that are now part of the standard amdlib package. The calibra- 
tion process includes an estimation of the calibrators' size and 
their uncertainties from various catalogs, a determination of their 
transfer function and their evolution during the whole night, and 
a computation of the calibrated visibilities and phases. The final 
errors on the measurements include uncertainties on the calibra- 
tors' diameter, the atmosphere transfer function fluctuations and 
intrinsic errors on the measurements. 



4. A qualitivative analysis of our dataset 

The Bry line for all the observed targets is clearly in emission 
(See Fig [2] and Table [3] for a summary of their main spectral 
characteristics). In almost all cases, the MR and HR data also 
exhibit a drop of visibility in the emission line caused by a vari- 
ation of the circumstellar environment extension and relative 
flux between the continuum and the line, as already explained 
in Meilland et al. (2007a). They also exhibit an "S" shape or 
more complex variations of the differential phase in the line, and 
the HR data clearly show that some visibility variations are "W" 



shaped (See e.g. Fig. 3). These characteristics are clear evidence 
of a rotating equatorial disk as described in Meilland et al. (201 1) 
for 6 Scorpii. The MR and HR data for all targets are presented 
inFigs.[3]to[T0] 

We note that for the HR observations executed before the re- 
placement of a disturbing optical element in front of the VLTI 
InfraRed Image Sensor (IRIS) in 2010, some instrumental mod- 
ulations with a very high frequency are seen in the HR data. This 
is especially the case for to Car. To enhance the data quality, 
we decided to filter these modulations with a standard Fourier- 
transform low-pass-filter technique. 

Some targets are also partly resolved in the K-band contin- 
uum. However, considering the uncertainties on the calibrated 
absolute visibilities, it is difficult to determine an accurate ex- 
tension of the circumstellar disk in the continuum for all targets. 
Therefore, we decided for efficiency to use only the differential 
visibility (i.e. visibility of each spectral channel divided by the 
mean visibility) for our kinematics study and the absolute con- 
tinuum visibility to determine the K-band extension when possi- 
ble. 

Among our sample we did not detect any new compan- 
ion, and 6 Cen remains the only star for which a companion 
was visible in the interferometric signal, as already described 
by Meilland et al. (2008). The individual comments for the ob- 
served objects are the following: 

• a Col : We obtained three measurements in HR and one in 
LR. The data S/N is quite high and the uncertainties on the 
differential quantities are on the order of a few percents. The 
target is clearly resolved in the line and the measurements in 
HR clearly exhibit the typical visibility and phase variations 
of a rotating disk with the major axis roughly perpendicu- 
lar to the polarization measurement obtained by Yudin 2001 
(see Table [TJ. The quasi-symmetric double-peaked line pro- 
file indicates that the object is seen from an intermediate in- 
clination angle and that no major inhomogeneity is present in 
the disk. The object is also partly resolved in the continuum, 
with a visibility of about 0.9±0.05 for the longest baselines. 
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Fig. 2. Bry spectra from our VLTI/ AMBER observations. The 
data were taken in HR mode for all observations expect for p 
Car and 6 Cen (MR mode). For k CMa the red and green lines 
represent the 2008 and 2010 observations, respectively. 

• k CMa : This star was observed for the first time in December 
2008 and we have obtained two measurements in HR and 
five in LR. It was observed again in January 2010 and two 
new measurements in HR were acquired. The S/N is quite 
high except for the last observations in 2010. All HR mea- 
surements at the two epochs show the typical visibility and 
phase variations of a rotating disk. The disk major axis seems 
to be roughly perpendicular to the polarization measure- 
ment. The double-peaked line profile and some visibility and 
phase variations are strongly asymmetric. These asymme- 
tries probably stem from an inhomogeneity in the disk that 
was already described in Meilland et al. (2007b). We note 
that the line profile violet-over-red (V/R) peaks ratio did not 
significantly vary between our 2008 and 2010 observations. 
The LR data show that the star is at least partly resolved in 
the continuum, but the calibration is not accurate enough to 
fully constrain the continuum extension. 

• to Car : We obtained two measurements in HR and one mea- 
surement in LR. Although the S/N is lower compared to 
a Col, a visibility and phase signals typical of a rotating disk 
are also present in the data. The major axis does not seem to 
be aligned or perpendicular to the polarization measurement. 
Because the line profile is double-peaked with a narrow shell 
line at its center, the star is probably seen from a high incli- 
nation angle. The LR data clearly indicate that the object is 
partly resolved in the continuum, i.e. V-0.85 for the longest 
baseline. We note that the second HR measurement cannot 
be calibrated accurately. 

• p Car : We obtained two measurements in MR, both show- 
ing a bright emission line, a drop of visibility, and "S" shaped 
phase variations in the line. Considering the phase amplitude 
and the lack of spectral resolution, it is hard to qualitatively 



Table 3. Bry line characteristics 



Star 


EW 


Peaks separation 


V/R 


Remarks 




A 


A (kms~') 






ffCol 


6.1 


9.9(138) 


1.04 


- 


k CMa(2008) 


16.8 


9.9 (138) 


1.13 




k CMa(2010) 


11.6 


10.4 (144) 


1.10 




a) Car 


3.5 


14.2(196) 


1.04 


Be-Shell 


p Car 


8.2 


-15(200) 


>1 


MR obs. 


6 Cen 


17.9 


-15(200) 


-1 


MR obs. 


yu Cen 


3.7 


10.4(144) 


1.03 




a Ara 


16.0 


15.1 (209) 


1.01 




o aqr 


12.6 


12.8(177) 


1.00 


Be-Shell 



determine the position of the major axis of the equatorial 
disk. We note that the profile and phase variations are asym- 
metric. This is probably cause by an inhomogeneity such as a 
one-armed oscillation, as proposed by Okazaki (1997). The 
object is barely resolved in the continuum (V-0.9 for the 
longest baselines). 

• £ Cen : Previous VLTI/ AMBER LR and MR observations 
published by Meilland et al. (2008) showed the binarity of 
this object. The authors found a contribution of the com- 
panion to the total K-band flux of 7% and a separation of 
68.7 mas. To check their results and constrain the system or- 
bit, the star was observed again in 2009, 2010, and 2011. In 
the 2009 MR data a bright Bry emission line, a drop of vis- 
ibility and "S" shaped phase variations are clearly visible. 
These data as well as the 2010 LR ones also contain modu- 
lations cause by the binarity of the object. We note that no 
obvious modulation is seen in the 201 1 LR dataset. This may 
mean that the separation is quite small, on the order of a few 
milli-arcseconds. 

• Cen : We obtained two measurements in HR. Because this 
target is quite faint, i.e., m^=4, the S/N is low. However, "S" 
shaped variations are still clearly present in the differential 
phases. We also detect a drop of visibility for one baseline, 
whereas it is clearly below the noise level for the other ones. 
The line profile is double-peaked and symmetric. The fact 
that the line intensity is quite low, i.e. 1.3 times the contin- 
uum, contributes to the weakness of the visibility drop. The 
object is unresolved in the continuum. 

• a Ara : This star was observed eleven times in LR mode in 
2007 but the data quality was too low to obtain more than 
an estimate of the disk extension (see Sect. 5.1). More re- 
cently, we have observed a Ara again during one full night 
and obtained seven measurements in HR mode. All data ex- 
hibit a clear signature of a rotating disk, with a major axis 
roughly perpendicular to the polarization measurement and 
compatible with Meilland et al. (2007a) results. The profile is 
double-peaked and symmetric. The object is partly resolved 
in the continuum for all baselines, but because of the unsta- 
ble weather conditions during the observing night, the un- 
certainties on the absolute visibility remains on the order of 
20%. 

• o Aqr : This is the faintest target of our sample with m^=4.6. 
We obtained one measurement in HR mode. The S/N is very 
low, and we were close to the instrumental sensibility limit 
considering the quite poor weather conditions during the ob- 
servations (seeing of 1" and a coherence time of less than 
3ms). However, an "S" phase signal is still clearly visible in 
the data, at least for the longest baseline. Clues of a small 
drop of visibility in the line are also present. The profile is 
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typical of Be stars seen at high inclination angle, i.e. double- 
peaked with a shell absorption line at its center. The object 
is clearly unresolved in the continuum. 

5. Modeling 

In this section we analyze our interferometric data using "toy 
models" of increasing complexity. In the first subsection the K- 
band continuum visibility is used to estimate the circumstellar 
environment extension. In the case of the Be binary star S Cen we 
tried to constrain the components' separation at various epochs. 
Finally, in the last subsection, the differential visibility and phase 
are used to constrain the circumstellar environment's kinematics. 

5.1. The disk extension in the continuum 

To estimate the circumstellar disk extension in the K-band con- 
tinuum with our interferometric measurements, we used a sim- 
ple two components model. Because all central stars are unre- 
solved or barely resolved even with the longest baselines (i.e. 
V>0.95), they were all modeled as point sources. The second 
contribution is the circumstellar environment. Owing to the large 
uncertainties on most of our absolute visibility measurements, 
we were unable to determine any flattening of the environment, 
and we decided to simply model this contribution as a circular 
Gaussian distribution. For the same reason, we decided to set the 
circumstellar environment relative flux (F env ) to the value deter- 
mined from the fit of the SED with a much higher precision than 
if it had been deduced from the interferometric data only (See 
Section 2 and Table 1). Finally, the two-dimensional intensity 
map describing our model is given by 

I(x,y) = (1 - F em ,)6(x,y) + I^expi-^f], (1) 

where x and y are the Cartesian coordinates, 5(x,y) is the Dirac 
function, and <x is the standard deviation of the Gaussian distri- 
bution. In the following, the Gaussian distribution is not defined 
by its standard deviation, but by its full width at half maximum 
defined by FWHM^2.35cr. 

For each dataset we separated the low-spectral resolution 
data from the medium or high-resolution data and fitted the 
above-described models. The results are summarized in Table[4] 
In most cases, we see that the accuracy of the LR data is higher 
than MR or HR, but that stronger biases affect them. Few cases 
were very problematic and no relevant fit to the data could be 
obtained. These special cases are marked as "n. c." for "not con- 
strained" in this table. 

Using the distance and stellar radius from Table [T| we 
were able to infer the extension in D*. These values, also pre- 
sented in Table [4] clearly show that the K-band continuum is 
quite compact. However, because of the large uncertainties on 
our measurements, we were unable to determine whether or 
not the K-band extension depends on any stellar parameters. 
Ne vertheles s, we determined a mean size of the environment, 
i.e. FWHM=2.2 ± 0.3 D* for the whole dataset. 

a Ara and k CMa were already observed with the 
VLTI/ AMBER and modeled by Meilland et al. (2007a, 2007b). 
For a Ara the authors found a mean continuum K-band ex- 
tension of 6.0±0.8 mas, which is significantly higher than our 
2.4+1.1 mas measurement. However, they used an uniform 
disk component for the environment. Because, the environ- 
ment is not fully resolved, we therefore have to apply a x0.87 



Table 4. The K-band continuum extensions of the disks. 



Name 


1 env 


Gaussian 


r W rllvl 


2 

Xr 




^11 ACQ ) 


mas 






a Col (HR) 


0.25 


1.3 ± 0.7 


1.5+0.8 


0.7 


a Col (LR) 


0.25 


1.0 + 0.2 


1.9+0.4 


2.0 


k CMa (HR) 


0.47 


n.c. 


n.c. 


n.c. 


K CMa (LK) 


r\ ah 
0.4 / 


1.0 ± 0.3 


3. /+1.1 


4.5 


to Car (HR) 


0.20 


n.c. 


n.c. 


n.c. 


to Car (LR) 


0.20 


1.7 + 0.5 


3.1+0.9 


1.9 


p Car (MR) 


0.45 


1.1 + 0.3 


2.0+0.5 


3.1 


fi Cen (HR) 


0.37 


n.c. 


n.c. 


n.c. 


a Ara (HR) 


0.56 


2.4+1.1 


3.8+1.7 


7.1 


a Ara (LR) 


0.56 


1.9+ 1.3 


3.0+2.1 


5.2 


o Aqr(HR) 


0.31 


n.c. 


n.c. 


n.c. 



correcting factor to their measurements to convert the uni- 
form disk diameter into a Gaussian FWHM. Nevertheless, this 
6.0x0.87=5.3±0.7mas is still significantly higher than our mea- 
surement (>2cr). In the case of k CMa the authors found that 
the mean K-band continuum extension was lower than 2.7 mas, 
which is compatible with our 1.0 + 0.3 mas extension. 

5.2. The binarity of 6 Cen 

We separately fitted the 6 Cen data with a uniform disk + a com- 
panion star model to compare it with Meilland et al. 2008. The 
two-dimensional intensity map is given by 



I(x,y) = F+6(x,y)+F comp S(x-Aa,y-A6)+^pn( ^ +y2 ),(2) 

nD env V D em> ' 

where F*, F comp , and F env are the stellar, companion, and envi- 
ronment fluxes, respectively, (Aa,A6) the components separation 
in Cartesian coordinates, and n(t) is the rectangle function de- 
fined by LI(t)=l for t<l/2 and II(t)=0 for t>l/2. 

The fit procedure is similar to that of Millour et al. (2009), 
i.e. a mix of Levenberg-Marquardt descent with a set of Monte 
Carlo initial parameters. The only free parameters are the com- 
panion coordinates. The other parameters were set to the values 
determined by Meilland et al. (2008): primary flux F*=0.41, en- 
velope flux F en v=0.52, companion flux F comp =0.07, and envelope 
diameter D env = 1 .6 mas. The results are presented in Table[5] 

We found that the position of the companion varied signifi- 
cantly between the four epochs separated roughly from one year 
each. It made an almost complete revolution around the main 
star during this period, pointing to a typical period of the sys- 
tem of about five years. To constrain the orbital elements sig- 
nificantly, the star should be observed again several times with 
long-baseline interferometry or speckles interferometry. These 
data should also be completed by radial velocity measurements. 

Nevertheless, we tried to determine a first estimate of a pos- 
sible orbit. It has a probably very low eccentricity, but a very 

Table 5. Evolution of the 5 Cen component separation. 
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PA 
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60.9 


-31.7 


68.7 


117.5 


2009-03 


-34.5 


0.79 


34.5 


-88.7 
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-73.8 


24.5 


77.8 


-71.6 
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Fig. 3. a Col differential visibility and phase from our four VLTI/ AMBER HR measurements (red line). Each row corresponds to 
one VLTI/ AMBER measurement (three different baselines). The visibility and phase of the best-fit kinematics model is overplotted 
in green. 
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Fig. 4. k CMa visibility and phase from our four VLTI/ AMBER HR measurements (red line). The two first measurements are from 
2008 (first two row) and the two other from 2010. The best-fit kinematics model is overplotted in green. 
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Fig. 5. a) Car visibility and phase from our two VLTI/ AMBER HR measurements (red line). The best-fit kinematics model is 
overplotted in green. 
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Fig. 6. p Car visibility and phase from our two VLTI/ AMBER MR measurements (red line). The best-fit kinematics model is 
overplotted in green. 
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Fig. 7. 6 Cen visibility and phase from our VLTI/ AMBER HR measurement (red line). The best-fit kinematics model is overplotted 
in green. 
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Fig. 8. fi Cen visibility and phase from our two VLTI/ AMBER HR measurements (red line). The best-fit kinematics model is 
overplotted in green. 
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Fig. 10. o Aqr visibility and phase from our VLTI/ AMBER HR measurement (red line). The best-fit kinematics model is overplotted 
in green. 
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Fig. 11. 6 Cen binary separation and possible orbit. The mea- 
sured separations are plotted as squares. The possible orbit is 
plotted as an orange solid line and the circles represents the mod- 
eled positions at the epochs of observations. 

large inclination angle. We managed to obtain a good fit of the 
orbit with the following parameters: semi-major axis of 80 mas, 
periastron in January 2008, period of 5.2 yr, no eccentricity, in- 
clination angle of 81°, u- 212°, and D=l 10°. We note that con- 
sidering the low number of measurements and the probable low 
eccentricity, this possible orbit, which we overplotted in FigfTTl 
may not be unique. 

5.3. The equatorial disk kinematics 
5.3.1 . A simple "toy" model 

To quantitatively constrain the velocity fields in the circumstel- 
lar environment of the observed Be stars using interferometric 
observations, we developed a simple two-dimensional kinematic 
model of a rotating and/or expanding equatorial disk. This model 
has already been used to model three classical Be stars (see 
Delaa et al. 2011 and Meilland et al. 2011) and one A[e] su- 
pergiant star (Millour et al. 201 1) and it is described in detail in 
Delaa et al. (2011). The model geometry is completely ad-hoc: 
the star is modeled as a uniform disk and the envelope emission 
in the continuum and the emission line as two elliptical Gaussian 
distributions of different FWHMs but with the same flattening 
due to a projection effect of the geometrically thin equatorial 
disk, i.e., / = l/cos(i), where i is the object inclination angle. 

The emission maps were then combined with a two- 
dimensional projected velocity map of a geometrically thin ex- 
panding and/or rotating equatorial disk. For each spectral chan- 
nel in the line, an iso-velocity map projected along the line of 
sight was then calculated and multiplied by the whole emission 
map in the line. Finally, the whole emission map for each wave- 
length consists of the weighted sum of the stellar map, the disk 
continuum map and the emission line map within the spectral 
channel under consideration (see Fig[12]for an example of emis- 
sion map obtained in a narrow spectral channel). 

The model parameters can be classified into four categories: 

1. The stellar parameters: stellar radius (/?*), distance (d), in- 
clination angle (/), and disk major-axis position angle (PA). 

2. The kinematic parameters: rotational velocity (V r ot) at the 
disk inner radius (i.e., photosphere), expansion velocity at 
the photosphere (VoX terminal velocity (Vco), and exponents 
of the expansion (y) and rotation (J3) velocity laws. 



3. The disk continuum parameters: disk FWHM in the contin- 
uum (a c ), disk continuum flux normalized by the total con- 
tinuum flux (F c ). 

4. The disk emission line parameters: disk FWHM in the line 
(«i) and line equivalent width (EW). 

The star distance is taken from van Leeuwen (2007) and F c and 

are derived from the fit of the SED (see Table Q]). The other 
nine parameters are free. 

If the disk is directly connected to the stellar surface, the ro- 
tational velocity (V rot ) should be equal to the stellar rotational 
velocity. However, in some cases, V rot may exceed the stellar 
velocity if the star is not critically rotating and some additional 
momentum is transferred to the circumstellar matter. Finally, we 
considered in our modeling that V 10t is free with a higher maxi- 
mum value equal to the critical velocity (V c ). 

For each target we computed several hundreds of models 
to constrain the parameters, determined the uncertainties and 
tried to detect any degeneracy or linked parameters. Owing 
to the large number of free-parameters, an automatic model- 
fitting method would have resulted in the computation of mil- 
lions of models. Moreover, we clearly know each parameter ef- 
fect on the visibility and phase variations (see Section 5.3.2.). 
Consequently, we decided to perform the fit manually. For all 
targets we could exclude models with significant expansion ve- 
locity of more than a few kms -1 . Consequently, we decided to 
set the expansion velocities to zero in all our models. We then 
tried to constrain the seven remaining parameters (i, PA, V rot , 
/J, a c , a\, and EW\). To reduce the number of computed models, 
we started with a qualitative estimation of the parameters from 
our interferometric data (especially for PA, i, a c , a; and EWi) 
and explored the parameter space with decreasing steps to con- 
verge to the^ 2 minimum. To check for other minima, we also ex- 
plored the full range of possible parameter space but with larger 
steps. Finally, the parameter values for the best-fit models are 
presented in Table [6] The corresponding differential visibilities 
and phases are overplotted in Figsl3lto[T0l 

The fit quality is very good for three targets observed in HR 
mode: a> Car, yu Cen, and a Ara, and good for the two stars ob- 
served in MR mode, i.e., p Car and 5 Cen. It is still satisfying for 
a Col (i.e. x%=4), although the visibility and phase of one of the 
baselines could not be fitted simultaneously with the other ones. 
In the case of k CMa, the fit is significantly worse (i.e.^ 2 =6.8). 
This is mainly because of the strong asymmetry of this object, 
which is not taken into account in our simple model. Finally, the 
data obtained on o Aqr seem to be insufficient to fully constrain 
the model for this object (i.e. <!)■ 

5.3.2. About the model parameters 

To check the consistency of our modeling, we tried to determine 
the effects of all model parameters on the visibility and phase 
variations through the emission line. Some of them strongly af- 
fect the interferometric observables and are consequently easily 
and unambiguously constrained with only a few measurements 
whereas others are more difficult to infer: 

• the major-axis position angle (PA) has a huge effect on the 
phase variation amplitude and the shape of the visibility 
drops, as already explained in Meilland et al. (2011). For a 
non-fully resolved disk, the amplitude of the "S" phase varia- 
tion is proportional to the baseline length, but it also strongly 
depends on its orientation. The amplitude is maximum for 
baselines aligned with the major axis and null for the one 
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Fig. 12. Illustration of our kinematic model. Left: two-dimensional project velocity map for a purely rotating disk with an inclination 
angle of 45°. The blue and red colors represent the positive and negative projected velocities, respectively. Center: global emission 
map in the Bry line. It is given by a simple elliptical Gaussian with an elongation cause by the projection onto the sky plane. Right: 
emission map in a narrow spectral channel centered around 2 1664 A with a resolution of 1.8 A. It is obtained using the velocity map 
and the emission map. See Delaa et al. (201 1) for more details. 



aligned with the minor axis. For baselines that overresolve 
the disk, the differential phase loses its simple "S" shape and 
secondary effects become visible (for example see the case 
of a Col in Fig. 4). The shape of the visibility variations also 
varies from "W" for baselines aligned with the major axis to 
"V" for baselines aligned with the minor axis. 

• the line equivalent width (EWi) is mainly set by the fit of the 
line profile. It needs to be corrected by taking into account 
the photospheric absorption line (see Delaa et al. 2011 for 
more details). 

• the disk FWHM in the emission line (ai) influences the drop 
of visibility amplitude along all baselines and the amplitude 
of the "S" shape variations. It can also be constrained know- 
ing that the phase variations lose this "S" shape for baselines 
that fully resolved the disk. This parameter also influences 
the double-peak separation: the larger the disk, the smaller 
the separation. 

• the disk FWHM in the continuum (a c ) is mainly derived from 
the measurements in the continuum presented in Table [4] 
However, it also indirectly influences the amplitude of the 
phase variations because it can modify the ratio between the 
coherent flux (for an unresolved object) and the incoherent 
one (for a resolved object). The phase variation is propor- 
tional to the photocenter shift only for an unresolved object. 
Therefore, the more resolved an object is in the continuum 
the smaller will be the phase variations. For example, in the 
case of a Ara, we did not manage to fit the differential phase 
with a disk in the continuum that extends to FWHM=3.0 mas 
as measured by the absolute visibility (see Table|4|, but only 
with disks smaller than 2 mas. This may be due to a trun- 
cation of the disk, as explained in Chesneau et al. (2005) 
and Meilland et al (2007a), making it depart from a simple 
Gaussian shape. 

• the inclination angle (/) has a significant influence on the vis- 
ibility drop amplitude for baselines close to the polar orienta- 
tion. It is well constrained by comparing equatorial and polar 
baseline measurements. It also has an effect on the double- 
peak separations since it influences the projected rotational 
velocity. 

• the rotational velocity (V rot ) mainly influences the double- 
peak separation and the line "width". The faster the disk ro- 
tates, the larger is the double-peak separation. 



• the exponent of the rotation law (J3) also influences the 
double-peak separation. With a higher value of /?, the ve- 
locity as a function of the distance drops quicker, and the 
double-peak separation is smaller. Therefore, it is quite hard 
to distinguish the effect of V rot and /3. However, f3 also in- 
fluences the shape of wings of the line. The higher the value 
of j6, the larger the line wings. But, as already mentioned in 
Delaa et al. (201 1), the wings of Be star emission lines can be 
highly affected by non-kinematic broadening cause by non- 
coherent scattering, so that it remains hard to set both V lot 
and f3 unambiguously. Nevertheless, value of f3 of less than 
0.3 yields variations that are too sharp, which are not realis- 
tic. 

6. Discussions 

6.1. The rotational rate of Be stars 

Using the inclination angle determined from our kinematic 
model and the v sin i and the estimate of the critical velocity V c 
from Fremat et al. (2005), we were able to determine the ro- 
tational velocity of the observed targets and constrain the ratio 
V/V c . These results are presented in Table [7] We found a mean 
ratio of V/V c =0.82+0.08. In Fremat et al. (2005), the rotational 
rate is defined by the ratio of the stellar angular velocity to its 
critical one : 

— = VReqc (3) 

&c V c Req 

where R eq and R eq c are the equatorial radii (in polar radii) for 
stars rotating at V and V c , respectively. Under the assumption 
of the Roche model, R eqc =1.5 and R eq =1.29 for V/V c ^0.82. 
Consequently, the mean rotational rate for our sample of Be stars 
is £2/Q c =0.95±0.02. This value is significantly higher than the 
one determined by Fremat et al. (2005), i.e. 0.88, from their fit of 
photospheric lines of Be stars taking into account gravity dark- 
ening effects as proposed by Townsend et al. (2004). However, 
our sample is much smaller than theirs and we note that the in- 
clination angle determined from our modeling agrees with the 
one determined by Fremat et al. (2005) within lcr for all targets 
except a Ara (1.3cr) and a Col (2<x). 

The values of Q/£2 C for each target are also presented in 
Table [7] We note that the upper value of the uncertainties for 
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Table 6. Best-fit parameters obtained from our axisymmetric kinematic model. 



Parameter 


unit 


a Col 


*CMa 


a) Car p Car 


(5Cen 


/j Cen 


a Ara 


o Aqr 


Global geometric parameters 


R* 


(Ra) 


5.8 


5.9 


6.8 6.0 


6.5 


5.5 


5.5 


4 


d 


(pc) 


80 


202 


104 148 


127 


155 


81 


133 


i 


(deg) 


35+5 


35+10 


65±10 70±10 


35+15 


25+5 


45+5 


70+20 


PA 


(deg) 


10 


25±10 


5±5 -25±10 


40+10 


80+15 


88+2 


120+20 


Global kinematic parameters 


V rot 


(kms- 1 ) 


350±10 


480±40 


300±20 400±30 


500+50 


510+20 


480+20 


400+50 


P 




0.5±0.1 


0.5±0.2 


0.45±0.1 0.45±0.1 i 


0.5+0.3 


0.5+0.1 


0.5+0.1 


0.5+0.2 


K-band continuum disk geometry 


F c 




0.25 


0.5 


0.2 0.45 


0.45 


0.37 


0.56 


0.31 




(D*) 


2.±0.5 


3.5±0.5 


3±1 2±0.5 


2+1 


<3 


<2 


<10 


Bry disk geometry 




(D*) 


5.5±0.3 


6.5±2 


6.5+1 11+2 


9+2 


4+1 


5.8+0.5 


14+1 


EWi 


(A) 


7.0±0.5 


13+2 


5.8+0.5 10+1 


19+2 


5.6+0.3 


14.5+1 


12+3 


xi 




4.0 


6.8 


1.1 2.5 


2.3 


1.3 


1.7 


0.8 



some stars are not physically reasonable since stars should 
not rotate above their critical velocity under the Roche model. 
However, we assume that the equatorial radius remains equal to 
1.5 for stars rotating above V c . In Fig. [13] the Be stars rotational 
rates are plotted as a function of the effective temperature to see 
if we find any correlation as proposed by Cranmer (2005). We 
did not detect any correlation between Q/D. c and T e g . This strik- 
ing effect could be real, or due to a small sample and our large 
uncertainties. 

Two stars in our samples have already been studied using the 
VLTI/ AMBER in MR mode: a Ara (Meilland et al. 2007a) and 
k CMa (Meilland et al. 2007b). The inclination angle found for 
a Ara, i.e. 45+5°, is roughly compatible with the previous esti- 
mate, i.e. 55° whereas, for k CMa the two estimates are clearly 
not compatible, i.e. 35±5° in this work and 60±10° in Meilland 
et al. (2007b). The main differences between these previous stud- 
ies and our current work is that they were conducted at the very 
beginning of the VLTI/ AMBER instrument. At this time, the un- 
certainties on the absolute visibility measurements were poorly 
known and had probably been underestimated by the authors. 
Unlike in the present work, these authors mainly used the abso- 
lute visibility measurements to determine the disk flattening and 
thus infer the inclination angle. Consequently, their conclusion 
that k CMa rotates at about half its critical velocity is probably 
biased. Nevertheless, with its strong asymmetry cause by the in- 
homogeneity in the disk it is still hard to accurately determine 
the k CMa rotational velocity. 



Table 7. Rotational rate of our Be stars 



Stal- 


v c 


V sini 


i 


V/Vc 






km s 


km s -1 


deg 






er Col 


355+23 


192+12 


35+5 


0.95+0.23 


99 +u M 


k CMa 


535+39 


244+17 


35+10 


0.80+0.31 


95 +0 16 


to Car 


320+17 


245+13 


65+10 


0.84+0.16 


96 +0 05 


p Car 


401+28 


285+20 


70+10 


0.76+0.15 


O Q9+0.06 

u - yz --o.ii 


6 Cen 


527+29 


263+14 


35+15 


0.87+0.41 


97+ 31 

-Q-32 


yu Cen 


508+32 


155+4 


25+5 


0.72+0.20 


90 + " m 
-fi-17. 


a Ara 


477+24 


305+15 


45+5 


0.90+0.17 


98 +!i *' 

-Q.Q5 


o Aqr 


391+27 


282+20 


70+20 


0.77+0.21 


93+ 06 



6.2. The equatorial disk extension 

All objects we studied are, at least, partly resolved in the Bry 
line. We managed to significantly constrain the extension of the 
line emission for all targets. We found Gaussian FWHMs that 
range between 4 to 14 stellar diameters. We found no correlation 
between the Bry emission and other properties or characteristics 
except for the double-peaked separation (see the next subsection 
for the discussion on the disk kinematics.). It seems that the size 
is independent of the stellar parameters and of the infrared ex- 
cess or line equivalent width. The mean FWHM of the Bry line 
emission is 6. 1 + 2.9 D* 

As already explained in Sect. 5.1, most of the targets are 
also partly resolved in the continuum with FWHM that range be- 
tween 1.5 and 3.7 D* and a mean FWHM of 2.2 + 0.3 D*. This is 
roughly compatible with K'-band CHARA interferometer mea- 
surements by Gies et al. (2007) on the Be stars y Cas (2.4 D*), <p 
Per (3.3 D*), f Tau (5.5 D*), and k Dra (4.3 D*). 

These typical K-band continuum and Bry line extensions are 
significantly smaller than the disk size measured in Ha. For ex- 
ample a multi-line spectro-interferometric study of the circum- 
stellar environment of the Be star 6 Sco published in Meilland 
et al. 2011 shows that the emission extension was 1.6 times 
more extended in H (9D*) than in Bry (5.5 D+) . Other nar- 
row band studies in Ha published in Tycner et al. (2004, 2006) 
and Quirrenbach et al. (1997) also concluded that the typical Ha 
extension were of the same order. 

As explained in Meilland et al. (2007a, 2007b) using the 
SIMECA code (Stee et al. 1994) and Tycner et al. (2007) using 
the BEDISK code (Sigut & Jones 2007), these interferometric 
measurements can be used to constrain the circumstellar envi- 
ronment physical parameters: mass loss, disk mass, and temper- 
ature and density distribution. In a forthcoming paper, we will 
use these numerical codes and the available R (CHARA/VEGA), 
K (VLTI/AMBER), and N (VLTI/MIDI) band measurements to 
draw a more complete picture of the Be star circumstellar envi- 
ronments. 

6.3. The equatorial disk kinematics 

For all our targets, the simple kinematic model reproduced our 
VLTI/AMBER measurements very well. Therefore, it is clear 
that most of the Bry line emission comes from the equatorial 
disk. The disk kinematics is dominated by rotation, with a rota- 
tional law close to Keplerian for all targets. The putative expan- 
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sion velocity is far below the detectability limit of the instrument 
(i.e. <10 kms -1 ). These results fully agree with previous kine- 
matics studies performed with the VLTI/AMBER (Meilland et 
al. 2007a, 201 1) or CHARA/VEGA (Delaa et al. 201 1). 

If the disks are in Keplerian motion and the stars are rotating 
significantly below their critical velocity, two questions remain: 

• What additional mechanisms provide the amount of energy 
needed to launch the matter from the stellar surface? 

• How does the ejected matter gain sufficient kinetic energy to 
accelerate up to the Keplerian velocity? 

In other terms, the matter needs to be accelerated both radi- 
ally and azimuthally. Lee et al. (1991) proposed that the disk 
could be formed by the effect of the gas viscosity that drifts 
the matter outward. However, in his theory, a source to supply 
angular momentum at the stellar surface is still needed if the 
star is not critically rotating. Non-radial pulsations as proposed 
by Osaki (1986) could be a good candidate. Recently, Cranmer 
(2009) proposed a theory in which resonant oscillations in the 
photosphere could inject enough angular momentum to spin up 
a Keplerian disk even for the slowest rotating Be stars (down 
to 60% of V c ). However, the question remains open and other 
mechanisms such as magnetism (Yudin et al. 2010), radiative 
pressure (Abbott 1979) or binarity could also contribute to the 
ejection of matter. To solve the question, the connecting layers 
between the stellar surface and the inner part of the equatorial 
disk should be carefully studied. This could be done by studying 
many absorption and emission lines formed close to the photo- 
sphere and by comparing their morphology. 

6.4. Departures from our simple model 

Despite the generally good fit of our interferometric data using 
the simple kinematic model, there are two significant departures 
from the model that need to be investigated: 

• The most important one is the case of k CMa that is clearly 
showing an asymmetric profile and visibility and phase vari- 
ations. These asymmetries can only originate from inho- 
mogeneities in the circumstellar environment. To determine 
whether or not these inhomogeneities can be modeled as 
one-armed oscillation as proposed by Okazaki ( 1 997) a dedi- 
cated model needs to be developed. Such a thorough analysis 
of k CMa data is out of the scope of this paper. We note that 
despite the lower resolution, the p Car data are also showing 
clues of the same kind of asymmetries. 

• The second most important departure from the model con- 
cerns the measurement obtained on a Col with the short 
polar baseline that cannot be fitted simultaneously with the 
other ones. This measurement shows that the environment 
is almost fully resolved with B^27 m for this orientation, 
whereas it is clearly less extended with a baseline closer 
to the equator. This could be a clue of a polar wind, as al- 
ready detected by Kervella & Domiciano de Souza (2006) 
on the classical Be star Achernar. However, a detailed study 
of the circumstellar environment of Be stars along their po- 
lar axis is needed to definitively answer the question of the 
detectability of polar winds in the K-band (see Stee 201 1 for 
a detailed discussion on the polar winds of Be stars). 

These two questions concerning the circumstellar environ- 
ments of Be stars will be studied in detail in some forthcoming 
dedicated papers. 
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Fig. 13. Rotational rate of our Be stars sample plotted as a func- 
tion of their effective temperature. Our measurements for the 
eight targets and their relative uncertainties are plotted as blue 
vertical lines. The orange box represents the mean rotational 
rate, i.e. 0.95±0.02. 

7. Conclusion 

In this spectro-interferometric survey of classical Be stars we 
managed to resolve all targets, constrain their extension in the 
Bry line and, for some of them, in the K-band continuum. Using 
a simple kinematic model of a purely rotating disk, we were 
able to successfully model all our data, showing that most of 
the Bry emission originates from the equatorial region. The disk 
is fully dominated by rotation, and the rotation law was found 
to be Keplerian or quasi-Keplerian for all targets. We were also 
able to significantly constrain the stellar rotational velocity with 
our estimates of the inclination angle. We found a mean rotation 
rate of Q/Q c =0.95±0.02. With such a high value of the rota- 
tional rate, only a small amount of energy is needed to expel 
matter from the stellar surface. This conclusion is quite different 
from the estimate by Fremat et al.(2005) obtained by modeling 
photospheric lines. 

We did not detect any correlation between the stellar param- 
eters and the disk properties. However, the uncertainties on the 
measurements remain high and our sample of Be stars needs to 
be extended to definitively answer the question on the physical 
process or processes responsible for the mass ejection and the 
dependence of the Be phenomenon on the stellar parameters. 

In a forthcoming paper, these data will be analyzed using the 
radiative transfer codes SIMECA and BEDISK to fully constrain 
the circumstellar environment. 
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